Exact solution for the energy density inside a one-dimensional non-static cavity with 

an arbitrary initial field state 
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We study the exact solution for the energy density of a real massless scalar field in a two- 
dimensional spacetime, inside a non-static cavity with an arbitrary initial field state, taking into 
account the Neumann and Dirichlet boundary conditions. This work generalizes the exact solution 
proposed by Cole and Schieve in the context of the Dirichlet boundary condition and vacuum as 
the initial state. We investigate diagonal states, examining the vacuum and thermal field as par- 
ticular cases. We also study non-diagonal initial field states, taking as examples the coherent and 
Schrodinger cat states. 
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I. INTRODUCTION 

The first works investigating the quantum problem 
of the radiation generated by moving mirrors in vac- 
uum where published in the 1970s decade (see Refs. 
BLUSH)- Moore 0], in the context of a real massless 
scalar field in a two dimensional spacetime, investigated 
the radiation generated in a cavity with a moving bound- 
ary. Imposing Dirichlet boundary condition to the field 
and also a prescribed law for the movement of a bound- 
ary, Moore obtained an exact formula for the expected 
value of the energy-momentum tensor, assuming the ini- 
tial field state as the vacuum. The field solution obtained 
by Moore is given in terms of a functional equation, usu- 
ally called Moore's equation (which was obtained inde- 
pendently by Vesnitskii ||), for which there is no general 
technique of analytical solution. Law [6| obtained an 
exact analytic solution for the Moore equation for a par- 
ticular resonant movement of the boundary (see also Ref. 
01)- Cole and Schieve proposed a numerical method to 
solve exactly the Moore equation for a general law of mo- 
tion of the boundary. Approximate analytical solutions 
of the Moore equation were also obtained, for instanc e, by 
Dodonov-Klimov-Nikonov 0] and Dalvit-Mazzitelli [lfj| |. 
With different approaches from those adopted by Moore 
and Fulling-Davies 0, perturbative methods where 
developed to solve the problem of a quantum field in the 
presence of a single moving boundary and also in 
oscillating cavities 

The first works studying the problem of particle cre- 
ation by moving mirrors with initial states different from 
vacuum were also published about thirty years ago 0], 
showing that the presence of real particles in the ini- 
tial state amplifies the phenomenon of particle creation. 
Considering a thermal bath as the "in" field state, the dy- 
namical Casimir effect has been investigated for the case 
of a single mirror [TIL EH EH > and also for an oscillating 
cavity 0, EH EH • The coherent state as initial field state 
has been considered for one moving boundary in Refs. 
(lH Ejj EE EH ! and the superposition of coherent sates 
in connection to the investigation of decoherence via the 



dynamical Casimir effect [20j] . Squeezed states have also 
been considered by three of the present authors [IH . 

Recently, several works have also investigated the in- 
fluence of different boundary conditions on the dynam- 
ical Casimir effect [U, El El HI- In this context, it 
was showed that, for the single mirror problem, Dirichlet 
and Neumann boundary conditions yield the same energy 
density radiated when the initial field state is symmetri- 
cal under time translations [l5|, EH ■ 

In the present work we investigate the time evolution 
of the energy density for a real massless scalar field in 
a two-dimensional spacetime, inside a non-static cavity, 
considering both the Dirichlet and Neumann boundary 
conditions. We extend to an arbitrary initial field state 
the exact method for calculating the energy density pro- 
posed by Cole and Schieve 0, HH in the context of the 
Dirichlet boundary condition and vacuum as the initial 
state. Considering an arbitrary initial field state, we 
show that the energy density in a given point of the space- 
time can be obtained by tracing back a sequence of null 
lines, connecting the value of the energy density at the 
given spacetime point to a certain known value of the 
energy at a point in the "static zone" , where the initial 
field modes are not affected by the perturbation caused 
by the boundary motion. We investigate diagonal states 
(for which the static energy density is invariant under 
time translation), examining in particular the vacuum 
and thermal field states. We also study non-diagonal ini- 
tial field states (in this case the static energy density is 
not invariant under time translation), investigating the 
coherent and the Schrodinger cat states. 

The paper is organized as follows. In Sec. [Til consider- 
ing Neumann and Dirichlet boundary conditions and also 
an arbitrary initial field state, we write the field solution 
and the energy density in the non-static cavity, general- 
izing the formula found in the literature 0, 0, which is 
valid for the Dirichlet boundary condition and vacuum as 
the initial field state. In Sec. IIIII we discuss the energy 
density in the static situation. In Sec. IIVI we obtain the 
energy density in the non-static situation given in terms 
of the static energy density. In Sec. [V]we summarize our 
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main results. 



II. ENERGY DENSITY: EXACT FORMULAS 

Let us start considering the field satisfying the Klein- 
Gordon equation (we assume throughout this paper h = 
c = ks = 1): (df — d%) tp (t, x) = 0, and obeying 
conditions imposed at the static boundary located at 
i = 0, and also at the moving boundary's position at 
x = L(t), where L(t) is a prescribed law for the mov- 
ing boundary with L(t < 0) = Lq, where Lq is the 
length of the cavity in the static situation. We con- 
sider four types of boundary conditions. The Dirichlct- 
Neumann (DN) boundary condition imposes Dirichlet 
condition at the static boundary, whereas the space 
derivative of the field taken in the instantaneously co- 
moving Lorentz frame vanishes at the moving boundary's 
position (Neumann condition): d x >ip'(t' ,x')\b oun dary — 
0. Using the appropriate Lorentz transformation, this 
boundary condition can be written in terms of quan- 
tities in the laboratory inertial frame as follows: 
d x >ij)'(t',x') = \(L(t)d t + d x )i>(t,x)\\ = 0. The 

L\ / J \x=L(t) 

Neumann-Neumann (NN) boundary condition imposes 
(t,x)) \ x=0 = and \{t{t)dt + d x ) V (t,x)] I = 0. 

L V / J I L(t) 

The Neumann-Dirichlet (ND) boundary condition im- 
poses [d x ip (t,x)]\ x=0 = and ip(t,L(t)) = 0, whereas 
the Dirichlet-Dirichlet (DD) boundary condition imposes 
tp (t, 0) = 0, and tp (t^LCt)) = 0. Considering the procedure 
adopted in Refs. [l|, 12], the field in the cavity can be 
obtained by exploiting the conformal invariance of the 
Klein-Gordon equation. The field operator, solution of 
the wave equation, is given by: 

oo 

$(t,x) =A[i+BV (0) (t,x)]+ Y, [a n ipn(t,x)+H.c.], 

n=l-2/3 

(1) 

where -0 (O) = [R{v) + R(u)} /2 0, and the field modes 
tp n (t,x) are given by: 



4> n (t,x) 



7^ } («)+7>r(«) , (2) 



with ip£\z) = e -i(«+«Tfl(2) ) u = t~ x, v = t + x, tmd 
R satisfying Moore's functional equation: 



R[t + L(t)] - R[t - L(t)] =2. 



(3) 



The operators A and B satisfy the commutation rules 



A,B 



B, a n 



= o m 



In Eqs. {TJ 



and ([2]) we introduce a notation which enables us to put 
into a single formula the solutions for the four boundary 
conditions considered in the present work. In this sense, 
for A = 7 = 1 and [3 = 0, Eqs. |TJ and © give the NN 
solution. The other three possible situations are recov- 
ered if we consider A = and: j3 = and 7 = i for the 



DD case; (3 = 1/2 and 7 = i for the DN case; /3=l/2 
and 7=1 for the ND case. 

Heareafter, considering the Heisenberg picture, we are 
interested in the averages (...) taken over any initial field 
state annihilated by B. In this context, we will write the 
exact formulas for the expected value of the energy den- 
sity operator T = (Tbo(i, x)). We can split T, writing: 



r = %. r 



^non-vac , 



where 



tt|7| 



71=1—20 



(n + (3) R ,2 {v)+R ,2 (u) 



and 



with 



^"non-vac — ^{a^a) + T(aa)i 



T( & &) = 92(v) + g2(u), 



(4) 



(5) 



(6) 



(7) 



(8) 



9i(z) 



7T ItT 



x Re 



n,n'=l-2/3 



[e< n - n ')* R W [i?'(z)] 2 (4a„-)}, (9) 



n,n'=l-2/3 



x Re | e -<»W+2/3)7rii( z ) [i? '( z )] 2 

>}■ (io) 

The term 7~ vac is the local energy density related to the 
vacuum state, which is divergent. Following Ref. 0], 
we adopt the point-splitting regularization method and 
obtain T vac (now redefined) as the renormalized local 
energy density: 



T vac = -/(«) - f(u), 



(11) 



where 



| 7 | 2 I R'" , 



3(/3 - f3 2 



R 



/2 



. (12) 

In the above equation the derivatives are taken with re- 
spect to the argument of the R function. 

We can also write the non-vacuum part of the energy 
density in the analogous notation: 



non-vac 



-g(v) - g(u), 



(13) 



3 



where g = —(31+52)- Note that T vac and T^ta) de- the spacetime variables. Despite the fact that T^ta) an d 
pend on I7I 2 , which has the same value for the considered 
boundary conditions. On the other hand 7/aa) depends 
on 7 2 , which can differ by a sign, depending on the situ- 
ation. For further analysis, it is useful to write: 



7^.^ can depend on time, from the principle of energy 
conservation the total energy is a constant in time for 
the static situation. In fact, for an arbitrary initial field 
state, the integration of Eqs. (TlT)|) and (f2T))) yields: 



where 



T = -h(v) - h(u), 



h = f + g. 



(14) 



(15) 



T (Sa) dx = UnMn ' 



T, ( A\ dx 

(aa) 



0, 



(22) 



(23) 



III. STATIC SITUATION 

Let us first examine the formulas in the preceding sec- 
tion for the static situation t < 0, when the Moore equa- 
tion is reduced to R{t + Lq) — R(t — Lq) = 2 . For this 
case the function R is given by [1( : 



R(z) = z/L . 



(16) 



The functions /, g, g\ and 32, now relabeled, respectively, 
as f( s \ g( s \ g± and g% , are given by: 



|7| 2 7T 

24L2 



2 - 303 -/3 2 ) 



(17) 



7« = 



<h 



2L 2 



x Re 



G?r+<?n (is) 



■ J2 y/(n + J) in' + J) 

n,n' = l-2/3 

{e<"-"> z/io (a,W)}, (19) 



2L 



n,n' = l-2/3 



X Re | e -<«+™'+2/3)7r Z /i (a„a n /)} . (20) 

Note that is a constant and T vac (in this case the 
Casimir energy density and relabeled as T cas ) is such 
that (see Refs. gHH): 



where ui n = 7r(n + (3)/Lq, and ftf n = (at a n ) is the initial 
mean number of particles in the nth mode. Then, for 
the static cavity, the function T^L can give contribution 
for the spacetime behavior of the energy density T (now 
relabeled as T^), but not for the total energy stored 

in the cavity. The total static energy — T^dx 
can be written as £^ = £ cas + X)« ^nNm where £ cas = 
Jq° Teas dx. Observe that, for a static cavity, for any 
initial field state considered, we have: £( S ) ND = £(s)DN 
and £(*)DD = £«NN 

For the initial field states such that the density matrix 
is diagonal in the Fock basis, so that (d n d n >) — and 
(d^an') =Af n 6 nn > we have: 

gi s) (z) = 0, d z [g[ s \z)] = 0, 0, [/»«(*)] - 0, fl « < 0, 

(24) 

leading to an uniform energy density in the static zone, 
which is invariant under time translation. For this case 
we have: 



y(a) DD _ j-{s 



NN 



T (s) DN = T {s) 



ND 



(25) 



This extends the equalities T = T C ^ S N and T C ^ S D = 
T C ^ S N , showed in Eqs. (JUJ), fr om the vacuum to any other 
diagonal field state. 

For initial field states with the density matrix having 
nonzero off-diagonal elements, the function g^ s ' can have 
spacetime dependence, what leads to an energy density 
non-invariant under time translation. For two cavities 
with different boundary conditions, but the same length 
Lq and values for (fit a„<) we have: 



T (s) DD _ --(a) NN T (s) DN _ T {s) ND 
i (ata) ~~ J (ata> ' 1 (aU) ~ i (ata) ' 



(26) 



<Y DD rr- NN 

cas cas 



24L 2 



<r- DN /y- ND 

' 1 cas cas 



48L 2 , 



, (21) 



where the superscripts DD, NN, DN and ND mean the 
types of boundary conditions considered in the calcula- 
tions. Note that, for mixed boundary conditions (in this 
case ND or DN), the Casimir energy is positive, origi- 
nating a repulsive Casimir force (see Ref. [26j|). On the 
other hand, g^ is, in general, spacetime dependent, what 
implies that T/^ta) an d %aa)^ m this static situation rc- 

(s) (s) 

labeled respectively as Xg ta \ an d are functions of 



On the other hand, for the same values of {a n a n >), we 
get: 

j-(s) DD _ _- 7 -(a) NN ^-(a) DN _ _sj-(s) ND ^-n 
(aa) (aa) ' (aa) (aa) ' v ' 

However, different and suitable choices for (a n a n i), for 
instance {d n d n '} 1 and (a„d„') 2 can produce: 



^-(s) DD _ ^(a) NN ^(a) DN _ j-(s) ND 



"(aa> 2 ' (aa) 1 

as showed in Sec. ITVBl Eq. 



(aa) 2 



(28) 



4 




FIG. 1: Boundary trajectories (solid lines). The dashed lines 
are null-lines separating the regions I from II, and II from III. 



IV. NON-STATIC SITUATION 

Now, let us examine the cavity in the non-static situ- 
ation (t > 0). The field modes in Eq. ([]} are formed by 
left and right-propagating parts. As causality requires, 
the field in the region I (v < Lq) (see Fig. [1]) is not 
affected by the boundary motion, so that, in this sense, 
this region is considered as a "static zone" . In the region 
II (v > Lq and u < Lq), the right-propagating parts of 
the field modes remain unaffected by the boundary mo- 
tion, so that the region II is also a static zone for these 
modes. On the other hand, the left-propagating parts in 
the region II are, in general, affected by the boundary 
movement. In the region III (u > Lq), both the left and 
right-propagating parts are affected. In summary, the 
functions corresponding to the left and right-propagating 
parts of the field modes are considered in the static zone 
if their argument z (z symbolizing v or u) is such that 
z < L a . We remark that Eqs. HU), dTTJ) , (fl9j) and (j20|) 
are valid not only in the static situation (t < 0) but also 
in the static zone (t ± x < Lq). 

For a certain spacetime point (t,x), the field operator 
ip(t,x) is known if its left and right-propagating parts, 
taken over, respectively, the null lines v = z\ and u = Z2 
(where z\ = t + x and z 2 = t — x), are known; or, in 
other words, ip(t, x) is known if R(v)\ v=Zl and R{u)\ u=Z2 
are known. Cole and Schieve Q proposed a recursive 
method to obtain the function R for a general law of mo- 
tion of the boundary, tracing back a sequence of null 
lines until a null line gets into the static zone where 
the R function is known. For a brief discussion of this 
method, let us assume that (t, x) belongs to the region 
III, and that the null line v = z\ intersects the moving 
mirror trajectory at the point [ti,L(ti)] (see Ref. [8|). 
We have R(v)\ v=Zl = R[h + L{t x )]. Using the Moore 
equation Q, we get R(v)\ v=Zl = i?(u)| u=tl _ L(tl) + 2, 
so that, as pointed by Cole and Schieve, "the value of 
R increases by 2 every time there is a reflection off the 
moving wall" [8|. If t\ — L{t{) < Lq, then the null line 



U = t\ — L(ti) is already in the static zone, so that 
we can write R(u)\ u = tl _ L ( t q = [*i - L(ti)]/L , and 
also R(v)\ v = Zl = [ti - L(ti)]/L + 2. On the other- 
hand, if t\ — L{t\) > Lq, we can draw another null 
line v = t\ — L(ti) intersecting the world line of the 
static boundary at the point [t\ — L(ti),0]. Then we 
get R(u)\ u=tl „ L{tl} = R(v)\ v=tl _ L{tl) , and can write 
R(v)\ v=Zl = i?(t>)|,j =fl _£( tl ) + 2. Assuming that the null 
line v = ti — L(t\) intersects the world line of the moving 
boundary at the point [£2, L(t 2 )], we get t2 + L(t 2 ) = ti — 
L(ti), and we can write: R{v)\ v=Zl = R(v)\ v=t2+L ^+2. 
In the same way, considering the Moore equation, we can 
write R{v)\ v=Zl = i?(u)| u=t2 _ L(t2 ) -I- 4, and so on, up 
to the null line considered is in the static zone, where 
the function R is known. In summary, we can write: 
R(z) = 2n + [z - 2J27=i L (U)}/L 0, where n is the 
number of reflections off the moving boundary, necessary 
to connect the null line v = z (or u = z) to a null line in 
the static zone. 

It is interesting to notice the behavior of the left and 
right-propagating parts of the field modes in the cavity. 
Taking into account the above analysis for the R behav- 
ior, for DD and NN cases (/? = 0) we get: ipn\v)\ v=Zl = 

(0) / % I (0) / \ (0) / \ I 

iph {U)\u=ti-L(ti), l?fi M |a=i 1 -£(t I )=W (») l^d-ifti) , 

and ( y 9l 0) (u)| t , = 4 1 _ L ( tl )=v3i 0) (ii)| u=t2 _ L(t2 ). We can see 

that the value of ip^ does not change under a reflec- 
tion off the static or moving boundaries. The relevant 
fact is the sequence of the spacetime points where the re- 
flections occur, because it will map an initial null line in a 
non-static zone, into a certain null line in the static zone, 
where the value of is determined. This last comment 
is also valid for mixed boundary conditions (/3 = 1/2), 

but, in this case, we observe that the value of ifn^ 2 ^ does 
not change under a reflection off the static boundary, but 
changes sign under a reflection off the moving boundary. 

The behavior of the energy density in a cavity, for the 
vacuum as the initial state, is described by the function 
/, which obeys the following equation proposed by Cole 
and Schieve 1231: 



/[< + £(*)] = f[t-L(t)]A(t)+B(t), (29) 



where 



A(t) 



l-L'{t) 
1 + L> (t) 



(30) 



Bit) = 



1 



L'» (t) 



127T [1 + L'(t)f[l-U(t)] 
1 L" 2 (t) V [t) 

"47T [1 + L'(i)] 4 [l-L'(t)] 2 ' 



(31) 



which enable us to obtain f(z), and therefore T va c, re- 
cursively in terms of its value / = /( s ) in the static zone 
(seeRefs. (HH). 

For an arbitrary initial field state, the behavior of the 
energy density in a cavity is described by the function 
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/ (related to the vacuum part) and a new function g 
(related to the non- vacuum part) proposed in the present 
paper. For the latter one, from Eqs. (|3|), j9|) and (fTTjf we 
obtain the following equation: 



g[t + L(t)]=g[t-L(t)] A(t), 



(32) 



which enable us to get g[z), and therefore T non -vac ; re- 
cursively in terms of its value g(z) = g^ s \z) in the static 
zone. 

Now, for a complete analysis involving vacuum and 
non- vacuum parts, let us write the equation for the func- 
tion h (Eq. I|l4])) as follows: 



h[t + L{t)] = h[t- L{t)]A(t) + B(t). 



(33) 



The procedure to find h(z), solving recursively the Eq. 
(f3"3")> . starts by setting z = t\ + L(t\) and tracing the 
null line v = t x + L(t x ). From Eq. ([33]), after the first 
reflection traced back on the moving boundary, we get 
the relation h(z) = h[t 1 — L (ti)] A(ti) + B(ti). Now, as 
explained above, we set t\ — L(t\) = t 2 + L(t 2 ), and, 
after the second reflection on the moving boundary, we 
get: 

h(z) = h[t 2 — L (t 2 )] A(t 2 )A(ti) + B(t 2 )A(ti) + B(ti). 

This process goes on, up to a null line gets into a static 
zone, where h = h 1 -^ — + g( s > is known. At the end, 
we get: 



where: 



h(z) = h^(z)A(z) + B(z) 

h^(z) = f^+g^[S(z)}, 
A(z) = U^AiU), 



(34) 



(35) 



(36) 



B(z) 



B{t n )W^A(U) 4 
+... + B(t 2 )A(h) 



Bit^n^AiU) 
\- 13(h). 



(37) 



The number n of reflections and the sequence of in- 
stants ti,...t n depend on z. In Eq. (|35p . we have 
z(z) = t n — L(t n ). From this process, we observe that 
the value of / does not change under a reflection off the 
static boundary, but it changes under a reflection off the 
moving boundary, as expected. We also note, from Eq. 
(l33|) . that the function h, after a reflection, changes in 
the same way for all boundary conditions considered in 
the present work (for the same boundary motion). We 
remark that, in Eq. (fM)) . the functions A and B depend 
only on the law of motion of the moving mirror, whereas 
all the dependence on the boundary conditions consid- 
ered in the present paper, and on the initial field state, 
are stored in the static region function h( s \z). We also 
observe that, for a generic law of motion, A and B are 
different functions, with the following properties: 

A{z) > V z, A{z < L ) = 1, B(z < L ) = 0, (38) 



which are directly obtained from Eqs. (|30|) and (|3Tj) . 

The energy density T can be given as in Eq. JJ), but 
now with the parts T vabC and 7^ n-vac rewritten as: 

T vac = -/ (s > \a(u) + A(v)] ~ B(u) - B(v), (39) 



gW[z(u)]A(u) -g (s) [z(v)}A(v 



(40) 



The instantaneous force acting on the moving bound- 
ary, in this two-dimensional model, is given by T(t 1 L(t)). 
From Eqs. (HH), and (gHJ) we see that two cavities 
with different boundary conditions and the same initial 
length presenting different vacuum and non- vacuum 
contributions for their initial energy densities, but with 
the same total energy density in the static zone, for 
a same (but arbitrary) law of motion will always exhibit 
identical evolutions in time for their energy densities T. 



A. Diagonal states 

Now, let us analyze and compare the behavior of 7^, ac , 
^non-vac and T for the initial field states for which the 
density matrix is diagonal in the Fock basis. From Eqs. 
(1241 and fiOl) we obtain: 



non-vac 



-fJ 



A{u) + A(v) 



(41) 



where g( s > is a constant factor. From Eqs. (f39|) and (|4Tj) 
we can see that all spacetime dependence in T comes 
from A and B. For this case we have: 



ij- DD _ rj- NN rj- DN _ <j ND 



(42) 



Since Eq. (|42[) is valid for any law of motion, it gen- 
eralizes the Eq. (|25|) . and because it is also valid for 
any diag onal state, it extends the conclusion obtained in 
Ref . [22| , where formulas correspondent to Eq. (|4"2"|) were 
obtained just considering the vacuum as the initial field 
state. 

The total energy stored in the cavity as a function 



of time can be written as £ (t) 
diagonal initial states, we have: 



^ {t) T(t,x) dx. For 



£(t) = -h^T 1 (t)-T 2 (t), 



(43) 



where Fi(t) = J^ (t) A{u) + A(v) dx and T 2 {t) = 



J*q B( u ) + &(v) dx. We emphasize that T\ and T 2 
depend only on the law of motion of the boundary, 
whereas all dependence on the initial field state or type 
of boundary condition is stored in the coefficient bS s >. 

Since the functions A and B are, in general, different 
one from the other, we conclude that 7^ ac , Tnon-vac and 
T exhibit, in general, different structures. On the other 
hand, one could ask about the conditions under which 



non-vac 



and T would exhibit the same structure 
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(the same curve apart from a positive scale factor or an 
additive constant) so that, for instance, if there are peaks 
and valleys they would be at the same positions in the 
cavity. One condition is provided by laws of motion for 
which B{z) and A(z) have a linear relation as given by: 



B(z) = kiA(z) + k 2 , 



(44) 



where k\ e fc 2 are constants. From the properties in Eq. 
\, we get k\ = — k 2 , resulting: 



(45) 



B(z) = h[A(z) - 1]. 
From ([21]), (@D} and 05]), we get: 

Tvac = -(J (S) + fcl)[i(«) + A(v)] + 2fc l' 



non-vac 



-gW[A(u) + A(v)], 



(46) 



(47) 



T 



(/(») + g W + ki)[A(u) + A(v)} + 2k x . (48) 



In the energy densities showed in Eqs. l[4"6" ]) -(|55 ]) , if 
the constant factors multiplying A{u) + A(v) are dif- 
ferent from zero and have the same sign, as the ratio 
a(z) = B(z)/[A(z) — 1] becomes more approximately 
equal to a constant value fci, more the structures of T vac , 
^non-vac and T become similar one to the other. Next, 
we exemplify these situations in the context of the vac- 
uum and thermal field, which are examples of a diagonal 
initial field states. 

For the thermal state we take into account that 
(aj,&„/) = 5 nn ,n{n,(3) and {a n a n >) = = 0, where 

n(n,(3) = l/(e' n+ ^' T — 1), with the Boltzmann con- 
stant (fcs) equal to 1. For this initial state, the energy 
density T non _ vac , now relabeled as Tt, is given by: 



7T|7l 

2Ll 



G(0,T) A(u)+A(v) 



(49) 



where Q(j3,T) = £n=i- 2/3 [(n + /3)/(e ( " +/3)/T - 1)] 
plays the role of a scale factor. Since (7(1/2, T) > G(0, T) 
(see Fig. [2]), for T > we have: 



T mixed ^_ t non-mixed 

It > It , 



(50) 



where the superscript "mixed" means ND or DN cases, 
whereas "non-mixed" means NN or DD cases. 

Note that if we set up two different cavities, one of 
them with non-mixed and the other one with mixed 
boundary conditions, both with a thermal field as the 
initial state with temperatures given, respectively, by T\ 
and T 2 such that 0(1/2, T 2 ) - Q{Q,T X ) = -1/16, we get 
for both cases the same initial configurations of energy 
densities in the static zone. Then for a same arbitrary 
law of motion for both cavities, they will exhibit the same 
time evolutions for the energy densities T (and hence the 
same force acting on the boundaries). An interesting case 
occurs if we set up T% » 0.34 and T 2 — 0. In this case we 




FIG. 2: Plot of Q{fi,T). The solid curve corresponds to 
£(1/2, T), whereas the dashed one to 5(0, T). 



obtain that the static and dynamical Casimir effect for 
the vacuum state in the mixed case (where the Casimir 
static force is repulsive) is mimicked by the non-mixed 
case with the mentioned specific temperature. 

Let us now investigate a particular movement de- 
scribed by: 



L (t) = L 



■ ( pnt\ 

£Sm UJ 



(51) 



where Lq = 1, p = 2 and e is a dimensionless parame- 
ter. This oscillatory boundary motion is investigated in 
several papers (see, for instance, Refs. d, |jj|)- In these 
papers the calculations were developed in the context of 
analytical approximate methods, considering small am- 
plitudes of oscillation (|e| <C 1). Here, we investigate the 
law of motion in Eq. (JST]) for e = 0.01 and e = 0.1. The 
former value is in better agreement with the condition 
|e| <C 1 than the latter one. Our intention now is to ver- 
ify, from an exact point of view, the similarity between 
the structures of 7^ a c, ^non-vac and T for both values 
of e. In Fig. [3] we plot the ratio a for the two values: 
e = 0.01 and e = 0.1. We observe that the ratio a is 
more approximately constant for e = 0.01 (dashed line) 
than for e = 0.1 (solid line), so that, since /W + fc 1 < 
and < 0, we expect that T VSiC , Tnon-vac and T exhibit 
more similarity for the former than for the latter value 
of e, for any diagonal initial state or any of the boundary 
conditions considered in the present paper. In fact, this 
is what can be visualized in Figs. 2] and El where we show 
7vac (dotted lines) and Tnon-vac (solid lines) for the law 
of motion (|5ip and the thermal bath with temperature 
T = 1 as the diagonal initial field state. In Fig. 2] for 
the case e = 0.01, both dotted and solid lines have their 
minimum and maximum values at the same positions, as 
predicted by the approximate analytical approach con- 
sidered by Andreata and Dodonov [IJI , which was based 
on e < 1. On the other hand, in Fig. [5j (case e = 0.1), if 
compared to the dotted line, the solid line has additional 
maximum points at x = 0.28 and x = 0.72, and also addi- 
tional minimum points at x = 0.42 and x = 0.58, so that 
differences between the structures of 7va C and 7~ non - vauC 
start to become more evident. 
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A A A A 




FIG. 3: The ratio a(z) = B(z)/[A(z)-l] for the law of motion 
given in Eq. (I51|l , The dashed line corresponds to e = 0.01, 
wheres the solid line corresponds to e = 0.1. 



FIG. 4: The energy densities 7Vac (solid line) and 7non-vac 
(dotted line) at the instant t = 1, for both DD and NN cases. 
The law of motion considered is showed in Eq. (|51[) , with the 
parameter e — 0.01. For the curve correspondent to 7non-vaci 
it is being considered a the thermal bath with temperature 
T = 1 as the initial field state. 



B. Non-diagonal states 

When we consider non-diagonal states, the function 
g("' in Eq. (|4"0|) can depend on z{z) or, in other words, it 
can depend on the null line in the static zone where it is 
being calculated. Since Eq. ([20)) depends on 7 2 , now T 
can be different for DD and NN cases (it is the same for 
DN and ND cases). We can see that T vac , T non . vac and 
T exhibit, in general, different structures. 

The total energy stored in the cavity as a function of 
time can be written as: 



£{t) = -f (s) T 1 (t)-T 2 (t)-T 3 (t), 



(52) 



where T 3 (t) = ^ {t) gW [z(u)]A(u) + gW [z(v)]A(v) 



dx. 

The function T 3 in Eq. (p)2"|) depends on the type of 
boundary condition, the initial field state and the law of 
motion considered. Next, we exemplify these situations 



FIG. 5: The energy densities 7Vac (solid line) and Tnon-vac 
(dotted line) at the instant t = 1, for both DD and NN cases. 
The law of motion considered is showed in Eq. (|51(l . with the 
parameter e = 0.1. For the curve correspondent to Tnon-vac, 
it is being considered a the thermal bath with temperature 
T — 1 as the initial field state. 



in the context of a coherent and a Schrodinger cat state, 
which are examples of non-diagonal initial field states. 

The coherent state is defined as an eigenstate of the an- 
nihilation operator: a n \a) — aS nno \a), where a = \a\e 
and n is related to the frequency of the excited mode 
[27l ] . The Schrodinger cat state [H, Hil is defined as a su- 
perposition of coherent states = N + e 1 ^ \— a)) , 
where |— a) has the same amplitude than \a) but with 
a phase shift of 7r, and the normalization constant N is 
given by N — (2 + 2e _2 l Q l cos</>) -1 / 2 . We combine into a 
single formula the results for the energy density T n on-vac 
(now relabeled as 7^) for both coherent and cat states, 
in the following manner: 



where 



T <ata) < 

1 (aa), 



9i 



(ata). 



(da) 



A{u) + A(v) 



(53) 



(54) 



<h 



(s) fz( U )]i( U )+ 5 ( s) [z(«)]i( W ), (55) 



with 
and 



2Ll 



(56) 



led 'Y^TT 

' ' 2 (no + (3) cos [2ir{n + [3)z/L - 26], 

(57) 

where C(0, \a\) = (1 — e _2 l"l" cos0)/(l + e~ 2 l"l coscj)). 
For 77 = we recover the coherent case, whereas for r\ = 1 
the cat state energy density is obtained. 

The two coherent states |a) and | — a) are not orthog- 
onal to each other and their overlap \ {a\ — a)\ 2 = e~ 4 '" 
decreases exponentially with \a\. For |a|=2 the overlap 
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is approximately zero [30(, then we expect that the cat 
state behaves like two coherent states. From Eq. ([53")) for 
\a\ « 2 we have C ~ 1 (see Fig. [6]), hence the energy den- 
sity for cat states and coherent states become the same. 
Both these energy densities are also the same, now for 
any value of \a\, if we consider the Yurke-Stoler state (a 
cat state with phase <fi = 

(s) 

Observe that g\ is constant in the static zone, whereas 
is spacetime dependent. From Eq. pTj) and l[5~4")) 
we see that Tnon-vac for diagonal states presents the 



same structure of 7/st 



An interesting result occurs 



when one considers the cat state with <fi = (odd cat 
state) and takes the limit |a| — > (leading to a sin- 
gle particle state) in Eq. (|53p . In this case we obtain 
gW = (n/2Ll) (n + (3) and T (aa) = 0. This means 
that the odd cat state in the limit |a| ^ presents its 
energy density with the same structure of a single particle 
diagonal state as expected. 

From Eq. ([57)) we see that 9 contributes to the shape 
of the energy density in the static zone. The Fig. [Jjshows 
the normalized coherent energy density T^\ n= o/\a\ 2 a s 
function of the position in a static situation, for fixed 9 
and several values of time and vice versa. Despite the 
fact that the energy density depends on 9, the 

total static energy £^ does not depend on this variable, 
as shown in Eq. (|2"5|) . On the other hand, since T^if) can 
depend on 9, from Eq. (f5"2"]) we see that the dynamical 
total energy £ can depend on this variable. 

For different boundary conditions, the same configura- 
tion of energy density can be obtained by choosing dif- 
ferent values of 9 in each case. Note that under a phase 
shift 9 — ► 9 + % we have the following symmetries: 



0) DD 0)NN, MDNi 0)ND, fKQ ^ 

9 2 \e=9 2 le+f > 9 2 \e = 9 2 le+f > ( 58 ) 




FIG. 6: C((j>, |q|) as function of |a|, for different values of i 




FIG. 7: The normalized energy density |^=o/|a| 2 (coher- 
ent state) with no = 1, as function of x/Lq. The solid, dotted, 
dashed, dashed-dotted and long-dashed lines can have differ- 
ent meanings, for instance: (a) In the DD case for t — 0, they 
represent, respectively, the energy density for 8 — 0, 7r/6, 7r/4, 
it/2 and 7r/3; (b) In the DD case for 9 — 0, they represent the 
time evolution of the energy density in the static situation 
(t < 0), corresponding, respectively, to t = 0, —1/6, —1/4, 
— 1/2 and —1/3. Also the correspondent NN situations can 
be obtained via Eq. (I59|l . 



thus we get: 



FINAL COMMENTS 



DD, 



(s) NN, 



C 



r/ s ) DN ( 2 )| 



( S ) ND, 



C 



9+f 
. (59) 

or, in other words, under the mentioned phase shift the 
energy densities for DN and ND (and also for DD and 
NN as shown in Fig. [Jj in the static zone are the same. 
Then, from Eqs. ([TJ| and we obtain that the energy 
density at any instant t satisfies: 



DD 



C 



it NN ^-DN 



it ND I 

= Tc \e+\ 



(60) 



In Fig. [8]we show the normalized coherent energy density 
7^/|a| 2 for DD, NN, ND and DN cases, obtained via 
numerical calculations based on the exact Eqs. (|5^|) and 
([55)) . In this figure, in which e = 0.01, the result for 
the DD case is in good agreement with the approximate 
analytical one obtained by Andreata and Dodonov for 
this boundary condition [1 91 ) . 



In summary, we extended to the case of an arbitrary 
initial field state the exact method proposed by Cole 
and Schieve [HI], and also took into account the Neu- 
mann boundary condition. The formulas obtained in the 
present paper enable us to get exact numerical results for 
the energy density in a non-static cavity for an arbitrary 
law of motion. Moreover, the structure of the formulas 
are meaningful. For instance, they exhibit that, apart 
from the law of motion, the evolution of the energy den- 
sity is completely determined by the total energy density 
in the static zone. In this sense we can set up two cavities 
with different boundary conditions and different initial 
field states, but in such manner that both have the same 
energy density in the static zone. In this situation we 
have (given the same law of motion for both cavities) the 
same time evolution for the energy density. We pointed 
that an interesting particular case occurs if we set up 
a non-mixed cavity with temperature 0.34 and a mixed 
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6000- 



5000- 



4000- 
p 3000- 



2000- 




0.26 



FIG. 8: The normalized coherent energy density 7^/|q| 2 for 
no = 1, t = 50.3 and the law of motion given in Eq. (|51[) 
with e = 0.01, as function of x/Lo- The solid line corresponds 
to both: T C DD with 9 = and T ( NN with 6 = tt/2. The 
dotted one corresponds to 7£ with 6 = and also T^ DD 
with 9 = it/2. The dashed-dotted line corresponds to Tq DN 
with 9 = and also T^ 13 with = tt/2. The dashed line 
shows T^ 13 with = and also T C DN with 9 = tt/2. 

cavity in the vacuum state. In this case we obtained that 
the static and dynamical Casimir effect for the vacuum 
state in a situation where the Casimir force is repulsive 
(the mixed case) are mimicked by another one for which 
the Casimir force is attractive (the non-mixed case) but 
with a certain non vanishing temperature. 

We obtained that the energy densities T vac , 7~ non - VSLC 
and T exhibit, in general, different structures. Wonder- 
ing about the conditions under which they would present 
the same behavior, we found that this particular situa- 
tion can occur for diagonal states and laws of motion for 
which the approximation given in Eq. (|45[) is valid. We 
showed that this condition is satisfied by the laws of mo- 
tion investigated in Ref. [l9| . where the same structure 
for these energy densities was obtained. On the other 
hand we showed that if the condition (|4"5)l is not satis- 
fied, these structures become different. 

Focusing on the influence of the boundary conditions, 
we obtained that, for a same diagonal initial state and 



(arbitrary) law of motion, the energy density for the DD 
case is identical to that obtained for the NN case, and the 
same occurs between ND and DN cases. This conclusion 
generalizes to any diagonal state the correspondent one 
obtained only for the vacuum case in Ref. [22(, and also 
extends to the cavity problem the conclusion that the 
force acting on the moving boundary is invariant under 
the change N D when the initial field state is invariant 
under time translations (in Refs. [HI, [l7| this conclusion 
was obtained for the one-boundary case). Investigating 
the thermal field, we obtained specific formulas and con- 
cluded that, for a same temperature of the field state in 
the static zone, the energy density for the mixed case 
is bigger than the one obtained for the non-mixed case. 
For non-diagonal states, which have the energy density in 
the static zone non-invariant under time translation, we 
obtained that the energy density for the DD case in the 
dynamical situation is in general different from the one 
obtained for the NN case, and the same is valid between 
DN and ND boundary conditions. However, we showed 
that in special situations the equalities can be recovered. 
Investigating the coherent and Schrodinger cat states, 
we obtained specific formulas for these cases and veri- 
fied that the difference between DD and NN cases (and 
also between ND and DN cases) is stored in T(aa) i part, 
whereas the T^d) ( P ar t behaves like the energy density 
of a diagonal state (all these conclusions are valid for any 
law of motion) . We applied our method of calculation to 
the specific situation of a coherent state in a DD cavity 
with an oscillatory motion with small amplitude, verify- 
ing that the approximate result found in the literature 
for the energy density [T^ | is in good agreement with the 
one obtained via the exact method presented here. In 
addition, we also showed the correspondent results for 
NN, ND and DN cases. 
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